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ABSTRACT 

Aims. We have studied the flare characteristics of 55 AGN at 8 dilferent frequency bands between 4.8 and 230 GHz. Our extensive 
database enables us to study the various observational properties of flares in these sources and compare our results with theoretical 
models. 

Methods. We visually extracted 159 individual flares from the flux density curves and calculated different parameters, such as the 
peak flux density and duration, in all the frequency bands. The selection of flares is based on the 22 and 37 GHz data from Metsahovi 
Radio Observatory and 90 and 230 GHz data from the SEST telescope. Additional lower frequency 4.8, 8, and 14.5 GHz data are from 
the University of Michigan Radio Observatory. We also calculated variability indices and compared them with earlier studies. 
Results. The observations seem to adhere well to the shock model, but there is still large scatter in the data. Especially the time 
delays between different frequency bands are difficult to study due to the incomplete sampling of the higher frequencies. The average 
duration of the flares is 2.5 years at 22 and 37 GHz, which shows that long-term monitoring is essential for understanding the typical 
behaviour in these sources. It also seems that the energy release in a flare is independent of the duration of the flare. 
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1. Introduction 

Active galactic nuclei (AGNs) exhibit variability across the 
whole electromagnetic spectrum, and it is often possible to see 
distinct flares in their flux density curves. By studying the dif- 
ferences in the variability of the various AGN types we gain un- 
derstanding into the physics involved in producing these events. 
The study of the radio regime is of special interest because it 
relates directly to the shock development in the jets. 

Variability studies of a large sample of AGNs at lower radio 
frequencies of 4.8-14.5 GHz have been done at the University 
of Michigan Radio Observatory (UMRAO), where AGNs have 
been monitored since the 1960s. Hughes et al. ( 1992) studied the 
long-term var iability time scale s of quasars and BL Lacertae ob- 
jects (BLOs). lAller et al.l(ll999l) studied th e total flux density and 
linear polarisation properties of BLOs and lAller et al.l (l2003h the 
variability of the complete Pearson-Readhead radio sample. 

At Metsahovi Radio Observatory AGNs have been moni- 
tored at 22 and 37 GHz for almost 30 years. Varia bility studies 
were first performed after 5 years of observations (Valtaoja et alj 
[1988; Valtaoja et al. 1992a,b) and updated by Lainelad99l 
with 12 yea rs of data. Long-term va riab ility time scal e s wer e 
studied by iLainela & Valtaoial (Il993h and iHovatta et all (l2007l) . 
Southern AGNs were studied at the hig her f requencies of 90 
and 23 GHz by [ Tornikoski et al.' (1993") and iTornikoski et al.l 
( 120001) . Also lOara mella et al. (2004) studied the variabihty 
properties of a sample of AGNs using the data from UMRAO 
and Metsahovi. 

Apart from these studies, most of the work done on the 
variability of AGNs concentrates on individ ual sources and 
their multifrequency behaviour. For example, in lPvatunina et al.l 



(l2006h and iPvatunina et al.l ( |2007|) seven AGNs are studied 
at 4.8-37 GHz. They decompose the flux density curves into 
Gaussian flares and study the properties of the flares and also 
the duty cycles in these sources. While this can describe and 
model one source very well, studies of large samples are impor- 
tant so that statistical diff'erences between the source types can 
be studied. 

The flares in the ra dio regime are explained by adia- 
batic shock-in-jet models dMarscher & Geadll985l: [Hughes et al.l 
1985). In these models an adiabatic shock is formed somewhere 
in the base of the jet and as it moves downward in the jet, we 
see synchrotron radiation at different frequencies and the ob- 
served properties depend on the phase of the shock de velopment. 
According to the model by iMarscher & Geaij (1 19851) the shock 
goes through three phases: The growth stage where Compton 
losses dominate and the flux density increases at all frequen- 
cies, a plateau stage where the energy gains and losses are equal 
and the peak frequency moves to lower frequencies, and a de- 
cay phase where adiabatic losses start to dominate and the flux 
density declines. 

I Valtaoja et al.l (Il992al) developed a parametrisation of the 
Marscher & Gear model - a generalised shock model - where 
the observed properties depend on the peak frequency of the syn- 
chrotron flare and how it relates to the observing frequency. If the 
peak frequency Vmax, where S{v^.^y,) reaches its highest value, is 
higher than the observing frequency Vobs, the flare is classified as 
a high-peaking flare. In this case the shock reaches its maximal 
development above the observing frequency and is in its decay- 
ing stage when observed. We should also see time delays towards 
lower frequencies. The flare is called low-peaking if Vobs > v^ax- 
Now the shock is in its growing stage when observed and flux 
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density curves of nearby frequencies trail each other closely and 
there should not be any time delays. 

This model can be tested by studying individual flares in 
AGNs. In this paper we used the Metsahovi monitoring sample 
as a starting point and extracted 159 individual flares from the 
flux density curves of 55 sources. We used 8 different frequency 
bands between 4.8 and 230 GHz to study the flare development. 
We also conducted a basic variability analysis by calculating the 
variability indices. 

The paper is organised as follows: In Sect.|2]we describe the 
source sample and data used. The methods used are described 
in Sect.|3] Section |4] includes the results of our analysis and the 
discussion follows in Sect.|5] The conclusions are drawn in Sect. 
|6] We use the value Hq - 72 km s~' Mpc~' throughout the paper. 
All the statistical analyses have been performed with the Unistat 
statistical package for Window^] (version 5.0). 

2. Data and the sample 

Our sam ple is select e d fro r n the Metsahovi mo ni toring 
source s dSalonen et al.1 119871; iTerasrantaetaTI 1 19921 1 19981 
12004 I2005h and from our monitoring sample of southern 
AGNs at the Swedish-ESO Submillimetre telescope (SEST) 
dTornikoski et al. 1996). The criterium for selection was that 
there must be at least one well-sampled outburst at two of the fre- 
quency bands between 22 and 230 GHz. This criterium was ful- 
filled for 55 sources. The minimum relative flux density change 
in these outbursts was 0.5 Jy at 37 GHz and between 0.4 and 
0.6 Jy at other frequency bands. In addition to the Metsahovi 
22, 37 and 90 GHz and SEST 90 and 230 GHz data we used the 
lower frequency data at 4.8, 8 and 14.5 GHz from the UMRAO 
database. Higher frequency data at 90, 150 and 230 GHz were 
also colle cted from the literature (Step pe et al.lI988llI992ilI993l: 
iReuter et al. 1997 ). Some of the southern sources are not moni- 
tored at Metsahovi and thus have data only at the higher 90 and 
230 GHz and the UMRAO low frequencies. 

Many of the sources have been monitored for over 25 years 
and here we used data from the beginning of our monitoring 
in the 1980s until 2005. Our study also includes unpublished 
data at 37 GHz from December 2001 to April 2005 and at 90 
and 230 GHz from the SEST from 1994.5-2003. The 37 GHz 
data for BL Lacertae objects from 2001-2005 are published in 
lNieppola"et al.l (l2007l) . The observation method and data reduc- 
tion pr ocess of Metsahovi data is described in iTerasranta et al.l 
(II998I) . 

Our sample consists of different types of AGNs, of which 14 
sources are BLOs, 4 are radio galaxies (GALs), 20 are Highly 
Polarised Quasars (HPQs), 12 are Low Polarisation Quasars 
(LPQs) and for 5 sources we had no information about their op- 
tical polarisation and they are therefore considered to be LPQs 
in this study . The BLOs are stud ied in more detail in a sep- 
arate paper (iNieppoIa et alj l2b08h . When studying the differ- 
ences between the source classes GALs are not considered be- 
cause of their smaller number in the study. The source sample 
is shown in Table [T] where the B1950-name, other name, clas- 
sification, possible EGRET detection, number of flares, number 
of datapoints since 1980 in each frequency band, maximum ob- 
served flux density at 37 GHz, variability index at 37 GHz and 
the median duration of flares at 37 GHz are tabulated. There 
is a clear selection effect in the sample as only sources with 
distinct flares were selected for the analysis. The radio prop- 
erties of these sources are interesting in view of the GLAST- 

' http://www.unistat.com/ 



satellite mission and therefore in Table [T] we have also marked 
the 29 sources that are high-confidence or possible detections of 
EGRET (.Hartman et al...l999.;.Mattox et al...2001,) . 

3. Methods 

We separated each well-monitored flare within the flux density 
curves and determined their start, peak and end epochs and flux 
densities for those epochs. Figures[T]and|9]show examples of flux 
density curves. In Fig. [T]the analysed flares of the LPQ source 
3C 273 (1226+023) ai-e shown at 37 and 4.8 GHz. Altogether 
159 flares fulfilled our criterium of being well-monitored in at 
least two frequency bands. From the start and end times we cal- 
culated the duration of the flares in the observer's frame, and 
using the peak epoch, also the rise and decay times at all fre- 
quency bands separately. In addition to the direct observational 
peak flux density of the flare, we used in our analysis the rela- 
tive peak flux density, which is the peak flux density minus the 
start flux density. Our flares are distributed between the classes 
so that there are 46 flares in BLOs, 15 in GALs, 59 in HPQs 
and 39 in LPQs. This means that on average every BLO has 3.3 
well-monitored flares during the 25 years while LPQs have one 
less (2.3). There are differences within the source classes as well. 
For example, the BLO source PKS 0735+178 has only one well- 
monitored flare during the 20 years' period it has been observed 
whereas BL Lac and OJ 287 both have 9 well-observed flares. 

We note that our definition of individual flares is based on 
visual inspection and is clearly different compared to models 
where a theoretical flare form is fitted to the observations (e.g. 
IValtaoia et all l999; Tu rler et alj|1999l) . Here we have, however, 
studied the observational properties of the flares and how they 
depend on each other. Therefore, a flare, as defined in this work, 
may include several individual shocks or events. 

4. Variability properties 

4. 1 . Variability indices 

We calculated the variability indices V for each source at all the 
frequency bands as defined in Eq.[T] 

^ _ (5 max - 0"5„„x) ~ min + O" s^^J 
(5 max - OS^J + (S min 

where 5 max is the maximum observed flux density of the source 
and (Ts^^ is its error and 5 min and os^.^ are the minimum flux 
density and its error The distribution of the indices in the dif- 
ferent source classes and at frequencies between 4.8 and 90 GHz 
are shown in Fig. |2] where also the median values are marked 
for each source class by vertical lines. 

We ran the Kruskal-Wallis analysis (KW-analysis) to see if 
there are any statistically significant differences in the variabil- 
ity indices of the frequency bands and source classes. The only 
difference between the frequency bands was that the 4.8 GHz 
indices differ from most other frequencies significantly within 
the 95% confidence limit. The indices at 4.8 GHz are smaller 
than those at higher frequencies. This can be understood as to- 
wards the lower frequency bands the flares are observed later 
due to synchrotron self-absorption. The flux densities are lower, 
timescales are longer, and it is also more difficult to distinguish 
the individual subsequent flares from each other. This is clearly 
seen in Fig. [T]if the 4.8 GHz flux density curve is compared to 
the 37 GHz curve. The relative flux densities are much lower and 
many individual flares seen in the 37 GHz curve, for example the 
flare in 1996, cannot be separated in the 4.8 GHz curve. 
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Fig. 2. Distribution of variability indices in different source classes and frequency bands. Median values for each source class are 
shown by vertical lines. 



When studying the distributions it seems that the BLOs have 
larger variability indices than the quasars. This is also confirmed 
by the KW-analysis which gives statistically significant differ- 
ences between BLOs and both quasar types at 14.5-37 GHz. At 
4.8 and 8 GHz LPQs and BLOs do not differ significantly but 
this is probably due to the low number of LPQs at these fre- 
quency bands. At 4.8 GHz there are only 6 LPQs and at 8 GHz 
only 5. These are also the more extreme LPQs with higher vari- 
ability. In general, the individual flares of relatively faint LPQs 
are often superposed at lower frequencies. 



4.2. Flare amplitudes 



The number of flares for which it was possible to determine the 
peak flux densities vary between the frequency bands. Most of 
the flares are defined at 22 and 37 GHz which were used for se- 
lecting the flares. We determined the observed peak flux density 
for 153 flares at 22 and 37 GHz whereas this was possible for 
only 33 flares at 230 GHz. The number of sources for which we 
were able to define the relative flux density is even smaller be- 
cause it was not always possible to determine when the flare 
started and what was the flux density then, even if the observed 
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Fig. 1. Upper panel: Flux density curve of the LPQ source 3C 
273 (1226+026) at 37 GHz. Flares used in our analysis are 
marked with arrows. Lower panel: Flux density curve at 4.8 GHz 

peak flux density is defined. In some cases we used the next min- 
imum after the flare for the relative flux density. 

The observed peak flux densities and relative flux densi- 
ties varied also significantly between the flares. The minimum 
peak flux density observed was 0.7 Jy and the maximum 56.9 
Jy. Clearly there are always some extreme objects which make 
the scatter large but the median values of these parameters seem 
to describe the overall behaviour in each source class and fre- 
quency band quite well. The median values of observed peak 
and relative flux densities of the flares are shown in Table |2] for 
each frequency band and source class. 

Considering the observed peak flux densities, there are no 
clear differences between the frequency bands according to the 
KW-analysis. When the relative peak flux densities are consid- 
ered the KW-analysis shows that only the 4.8 GHz dififers from 
the other frequency bands. Also the higher frequency bands from 
14.5 GHz to 230 GHz form one group while the 8-22 GHz fre- 
quency bands form another group. The median values have an in- 
creasing trend in the relative peak flux densities as the frequency 
increases. It also appears that in BLOs and HPQs the flux den- 
sities are somewhat higher at higher frequencies but the differ- 
ence is not large. For LPQs, the peak flux densities at the lower 
4.8-14.5 GHz frequencies clearly differ from those at the higher 
frequency bands but this is probably due to the much smaller 
number of sources at these frequencies. As was already noted in 
the previous section, at lower frequencies the flares are more su- 



perposed and it is the extreme objects that stand out. The fainter 
sources do not show distinguishable flares at lower frequencies 
and therefore are not included in the calculations. This is also 
seen in the KW-analysis which shows that at 4.8 and 8 GHz the 
LPQs differ from the other classes significantly. At 14.5 GHz 
the difference with the higher frequencies is not as clear when 
the relative peak ffux densities are considered. At the higher fre- 
quencies there are no significant differences between the source 
classes. 

4.3. Flare intensity and duration 

In some cases we could not calculate the duration of the flare 
because either the starting or ending point was not determined. 
This varied between the frequency bands so that we were able to 
calculate the duration for 150 flares at 37 GHz and for 31 flares 
at 230 GHz. The durations of the flares varied from 0.3 years to 
13.2 years. 

The duration seems to decrease with increasing frequency 
but again the difference is not very large, only 0.6 years at most 
between 4.8 and 90 GHz. According to the KW-analysis there 
are also no significant differences between the source classes 
even though the LPQs seem to have slightly longer durations. 
This could indicate that there are differences in the physical con- 
ditions of the jets in LPQs and other classes. Of particular note 
is the long median duration of the flares: They last about 2.5 
years at 22 and 37 GHz. This means that short multi-wavelength 
campaigns are not able to reveal the complete behaviour of these 
sources at radio frequencies. Long-term monitoring is needed to 
study the variability at these frequencies. 

The rise and decay times appear to be approximately equal 
and so the flares are relatively symmetric. At 22 and 37 GHz 
the median rise times are 1 year and decay times 1.3 years in 
accordance with the results of Valtaoja et al. ( 1999). There are 
no significant differences between the frequency bands but the 
decay times of the LPQs seem to be slightly longer which is also 
seen in the KW-analysis. 

In Fig. |3] we plotted the observed peak flux density of 
the flares against the duration for all sources at 37 GHz. The 
Spearman rank correlation for these parameters is r = 0.15 
with probability p - 0.03. Thus the correlation is not strong. 
At lower frequencies the correlation was slightly higher, de- 
creasing towards higher frequencies. If only quasars (HPQs and 
LPQs) are considered, there is a correlation of r = 0.29 with 
p = 0.002 at 37 GHz. The result is similar if HPQs and LPQs 
are tr eated separately. For B LOs there are no significant correla- 
tions (iNieppola et alJl2008h . 

There is a slightly higher positive correlation if the rela- 
tive peak flux density is studied against the duration but the 
distribution of points looks very similar At 37 GHz r - 0.26 
with p = 0.001. Again the correlations at lower frequencies 
are slightly higher, decreasing towards higher frequencies. This 
time the correlation of quasars is even higher with r = 0.36 and 
p - 0.0002. Also in this case there are no correl ations between 
these parameters if only BLOs are considered dNieppola et al.l 
2008). 

There is a selection effect which could affect the results. 
There are fewer weak flares of long duration. If such a flare is 
seen on the flux density curve, it may get overlooked as it re- 
sembles more a secular baseline flux density trend than a flare. 
Also the shorter the flares are, the more they show up like spikes 
in the flux density curve and are taken as flares more easily. 

One other significant effect is Doppler boosting. When the 
sources are more boosted, their flux density levels are higher 
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Table 2. Median values of duration, observed peak flux density and relative peak flux density for all frequencies and source classes. 
Number of flares used to calculate the value is also shown. 
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where t is the duration. So is the flux density, Dl is the lumi- 
nosity distance, z is the redshift and Dyar is the Doppler boost- 
ing factor taken from Hovatta et al. (in preparation) where also 
other jet parameters such as jet speed are studied in more de- 
tail. When the corrected durations are studied there are still no 
clear differences between the frequency bands according to the 
KW-analysis. The LPQs have slightly longer flares but the differ- 
ence is smaller than in the uncorrected case. Instead, at the lower 
4.8-14.5 GHz frequencies the BLOs differ from the other source 
classes due to their shorter flares. The corrected peak luminosity 
against the corrected duration is shown in Fig. |4] The positive 
correlation disappears and there are no significant correlations 
according to the Spearman rank correlation. 

4.4. Correspondence to the generalised shock model 

Accord ing to the generalised shock model (IValtaoja et al.l 
Il992al) . the observed properties of a flare depend on the ob- 
serving frequency and whether the peak frequency of the flare 
is higher or lower than the observing frequency. We studied how 
well these 159 flares correspond to the model by calculating a 
few parameters for each flare. We point out, however, that in this 
work we have not extracted separate shock events from flares but 
treated each longer outburst event as an individual flare. 

First we determined the frequency band in which the flare 
occurs first and we have given it rank number 1. The second 
frequency band then gets number 2 and so on. For example, 
if the flare peak is first at 90 GHz, then at 22 GHz and then 
at 37 GHz, the rank numbers are 1 for 90 GHz, 2 for 22 GHz 
and 3 for 37 GHz. In this procedure we ignored flares which are 
clearly multi-peaked and for which it is impossible to determine 
the actual peak at all frequency bands. We have not ignored fre- 
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Fig. 3. Observed peak flux density against the duration of the 
flare at 37 GHz. 



and at the same time the flares appear shorter Therefore we 
corrected the peak flux density and dur ation for redshift and 
Dopp ler boosting using Eqs.l2land[3](e.g. lKembhavi & Narlikarl 
119991) 
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Fig. 4. Corrected peak luminosity at 37 GHz against the duration 
of the flare corrected for redshift and Doppler boosting. 
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Fig. 6. Relative flux density rank number plotted against the fre- 
quency. Frequency with highest relative flux density has been 
ranked 1. 
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Fig. 5. Peak time rank number plotted against the frequency. The 
frequency band having the flare peak first has been ranked 1 . 



quencies with larger gaps and therefore we expect scatter in the 
results. We plotted the peak time rank numbers against the fre- 
quency in Fig. |5] in a bubble plot. The bigger the bubble, the 
more cases there are in that rank number and frequency. The 
plot shows that, as expected, the flares are peaking first at higher 
frequencies but there are also outliers and scatter. 

Another prediction of the model is that relative flux density 
should be highest at the peak frequency, which in the case of 
high-peaking flares should be in the upper part of our frequency 
range. We made a similar bubble plot of the relative flux den- 
sity rank number against the frequency (Fig. |6]l. In this case we 
have given rank number 1 to the frequency at which the relative 
flux density is highest. Rank number 2 is given to the frequency 
which has the second highest relative flux density and so on. Fig. 
|6]shows that the highest flux density is usually at 22 and 37 GHz, 
although the number of flares at higher frequencies is lower and 
therefore the result could be slightly biased. Also the sampling 
of the higher frequencies is more sparse and therefore the actual 
flare peaks could have been missed in the observations. 



A more straightforward way of studying the correspondence 
between the model and the observations is to normalise the rel- 
ative flux density to the flux density of the peak frequency Vmax 
of the flare. In our analysis we have simply taken v^^^ as the fre- 
quency band at which the relative peak flux density is highest. A 
more accurate way would be to interpolate the peak frequency 
but we find it unsuitable for our data which has large gaps es- 
pecially at the higher 90 and 230 GHz frequency bands, making 
the estimation of the true peak time and peak flux density of the 
flare difficult. According to Valtaoia et al. ( 1992a) the maximum 
flux density of the flare should correspond to the plateau stage in 
the generalised model. Frequencies higher than v^ax correspond 
to the growth stage and should have A^max v""»". Frequencies 
lower than Vjnax should have AS'max °^ where b dep ends on 
the mod el. This appr oach wa s originally used by Valtaoj aet al.l 
(Il992ah and later by iLaineial d 19941) . They used a smaller sam- 
ple of sources and as the sources had been monitored for only a 
short time, they used the overall minimum and maximum values 
ever observed as the relative flux density. Our approach should 
be more accurate as the flares were individually separated and 
the number of flares is much higher. 

In Fig. Q we plotted the normalised relative flux density 
against the normalised frequency. All the values were nor- 
malised to the frequency at which the relative flux density is the 
highest. We also fitted straight lines to below and above the peak 
frequency. We can see that there is large scatter in the values but 
th e trend is very much a s expected when compared to the model 
in Valta oia et al.l (Il992ab . The fit gives A^^ax v"**' for the ris- 
ing part and AS max v^^-^'^ for the declining part. We also fitted 
the data separately for quasars and HPQs and LPQs and the re- 
sults were very similar to those of the whole sample. According 
to the model, there should be a plateau between the slopes which 
is not seen in our analysis. We also tried to fit a plateau stage into 
the model, but two linear components describe the data more ac- 
curately . The same result was already obtained in lValtaoja et al.l 
(Il992al) . They got a slope of 0.52 for the rising part and -0.26 
for the declining part, which are very close to the values ob- 
tained here. iStevens et a l. ( 1994) also obtained a value of -0.28 
for the dechning part in their study of 17 blazars between 22 and 
375 GHz. 
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Fig. 7. Observed maximum relative flux density scaled to the fre- 
quency at which the maximum occurs. The two straight lines 
give best fits to the section below and above the peak frequency. 
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frequency against the normalised frequency. 



Time delays between the frequencies are also important but 
more difficult to study. This is because the sampling is often 
sparse and the actual peak time could have been missed during 
observations. If a source is observed once in two weeks, which 
is quite a common sampling rate for these sources at Metsahovi, 
the peak can occur at any time within a month. Often the gaps are 
longer due to weather constrains or maintenance. Therefore any 
time delays shorter than a month are not reliable. At the higher 
90 and 230 GHz frequencies the situation is even worse because 
the flares are faster and the data acquisition is even more sparse. 

We calculated median time delays at diff'erent frequency 
bands and for 37 GHz it is 34 days, increasing towards lower fre- 
quencies and being 181 days at 4.8 GHz. According to the gen- 
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Fig. 9. a) Flux density curve of the BLO source 0235+164 at 
37 GHz with each flare analysed marked with arrows, b) Flare 
number 3 in more detail, c) Spectra of every epoch marked in 
Figb). 



eralised shock model, there should not be significant time delays 
above the peak frequency but below it they should increase. We 
calculated the time delays between the peak frequency and other 
frequencies and plotted them against the normalised frequency 
in Fig. [8] It is clear that the trend is correct with time delays be- 
ing larger at frequencies lower than the peak frequency. There 
is also large scatter in the values which is expected because the 
time delays are so difficult to determine especially at higher fre- 
quencies. In this analysis we took into account all flares where 
the peak was not clearly dififerent at diff'erent frequencies and 
we did not ignore flares with poor sampling. Even time delays 
as large as 300 days can be due to p oor sampling at 230 G Hz. 
However, if compared to the model in Valtaoia et al.l (Il992a l). we 
see clear correspondence between the model and observations. 

We also studied one individual source, 0235+164, and its 
flare in 1993 in more detail. Figure |9^, shows the flux density 
curve of the source at 37 GHz, with every flare analysed for 
this source marked using arrows. Fig.|9j) shows the 1993 flare 
in more detail at 37 GHz. We have marked seven epochs of the 
flare and constructed simultaneous spectra at every epoch at all 
the frequencies used in the analysis (Fig. |9};). The spectra are 
simultaneous within two weeks but the 230 GHz points are in 
general more uncertain because of the more sparse sampling. 
The flare follows the shock model quite well with the flux den- 
sity first rising at all frequencies. When the flux density starts to 
decay the peak of the spectrum moves to lower frequencies. No 
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clear plateau is seen in this flare which agrees with the general 
results for the whole sample. 

5. Discussion 

We have studied the observational properties of radio flares for 
a large sample of densely monitored variable AGNs. For many 
of the sources the data sets span more than 25 years. For this 
study we only chose objects with clearly distinguishable radio 
flares, which biases this sample towards sources with notable 
variability and flares with relatively simple morphology. 

Within our sample, the BLOs have higher variability indices 
than quasars. This is in part caused by the fact that only the BLOs 
with truly distinct variability features were includ ed in this study. 
Compared to the variability indices obtained in Nieppo la et al.l 
(2007), our indices are higher, the media n variability in dices o f 
BLOs being 0.78 in fliis study and . 33 inl^ieppola et alj(l2007l) . 

The sample of iNieppola et al.l (12007 ) also includes fainter 
BLOs for which no well-structured radio flares could be ex- 
tracted from the flux density curves. The differences between 
the variability behaviour of these two BLO samples can be ex- 
plained by the fact that only a small minority of the full BLO 
sample actually shows any extreme behaviour in the radio do- 
main, most of them being faint and relatively non-variable even 
over several years of monitoring observation s. However, it is 
noteworthy that the large BLO sample of Nieppola et al.l (l2007h 
only includes 5 years of data compared to our 25 years, and at 
least some of the sources may radically change their behaviour 
over longer timescales, as discussed in Hovatta et aL (20Q7). 

Earlier studies of higher 90 and 23 GH z variability in- 
dices w ere done bv lTornikoski et aLl(ll993l) and iTornikoski et al.l 
( I2OOOI) . Their sample included southern AGNs observed with the 
SEST telescope. The variability indices are significantly lower 
than what we have obtained (for comparison we calculated the 
variability index defined as (^max - 'Smin)/'Smin for our sample) 
but also the number of datapoints is much smaller. The aver- 
age num ber of datapoint s in ou r sample is almost five times 
larger. In'Torniko ski et al] (l2000l) the effect of the number of dat- 
apoints was studied and they noticed that the variability indices 
increased with increasing number of datapoints, confirming our 
findings where only very long-term monitoring reveals the true 
minimum and maximum flux densities of most sources. Even 
though the highest flux densities during flares are typically ob- 
served in the 37 and 22 GHz frequency bands (as seen in Fig.|6]), 
the 90 GHz variability indices obtained in this paper are of the 
same order as those in the lower frequency bands, indicating that 
strong flares can indeed occur also in the mm-domain. 

When we studied the duration and intensity of the individual 
flares, it was evident that most flares are very long lasting. The 
median duration of flares at 22 and 37 GHz is 2.5 years, which 
means that short-term observing campaigns are rarely helpful 
for describing the "typical behavio ur" of these s ources. When 
we add to this the result from lHovatta et alj (l2007h that for these 
sources flares occur on average every 4-6 years, we can con- 
clude that monitoring campaigns lasting at least ca. 5-7 years 
are needed in order to catch a typical source both in a quies- 
cent and flaring state and to follow the flare evolution from the 
beginning to the end. 

The flare intensities increase with flare duration only to some 
degree. For flares lasting 2-3 years we can see almost as high 
peak flux densities (observed peak flux density as well as the rel- 
ative flux density during the flare) as for those lasting a few years 
longer. A striking result is that when we calculate the Doppler- 
corrected peak luminosities of these flares, the weak positive cor- 



relation disappears altogether and there is no correlation between 
the Doppler corrected peak luminosity against the Doppler cor- 
rected flare duration. This means that the energy release in a flare 
does not increase with duration for these flares. 

When studying the rise and fall times of the flares we see that 
flares in general have semi-symmetric shapes, with the decay 
time typically being 1.3 times the rise time. The morphologies 
of flares are typically not very simple, though, and flares with 
multiple peaks or other finer structure complicate the analysis. 
In general, flares at 37 GHz or lower frequencies rarely consist 
of smooth rise and decay components only, and a more thorough 
analysis is needed to fully understand the physical characteristics 
of radio flares and to separate the individual flare components 
(Hovatta et al. in preparation). 

Based on the observational results of this study, we have 
studied the correspondence between the observations and the 
generalised shock model bv lValtaoia et al.l (Il992al) . For 99 flares 
we had an estimate of the time delays between at least 3 fre- 
quency bands. Nearly half of them (46 flares) seemed to be- 
have as the shock model predicts with no time delays above 
the peak frequency and increasing time delays below it. In ad- 
dition in 34 cases the time delays were shorter than the sampling 
rate and therefore we cannot say whether they follow the shock 
model or not. Only in 20% (19 cases) the time delays and ampli- 
tudes did not follow the shock model. In these cases, however, 
the flares were defined to have multiple peaks and therefore the 
time delays were difficult to determine. Also the sampling rate 
may have been poor In any case the flares seem to follow the 
shock model in general, even though there is large scatter due to 
poor sampling and complicated structure of the flares at different 
frequency bands. Of the well-determined flares where we could 
study if they are high- or low-peaking, 25 were high-peaking 
and 3 low-peaking at 22 and 37 GHz. This indicates that the 
peak frequency of the flares in these sources is some where in 
the millimetre domain. This is in agreement with IStevens et al.l 
il994^ who found that in their sample of 17 blazars studied at 
22-375 GHz the peak was at 90 GHz or below. Our results for 
the whole sample are very similar to those for BLOs studied in 
Nieppola et al. (2008). 

When we studied the shape of the flare spectra normalised 
to the peak frequency, we found the best fit to be AS max v" "*' 



for the rising part and A5max 



for the declining part. 



We did not find a plateau between the slopes which is contra- 
dictory to the model but in agreement with ea rlier studies by 
iValtaoia et all (Il992ah and lLindfors et all ( l2006h . The decHning 
slope a - -0.24 corresponds to the optically thin spectral index 
ffthin for a pow er-law electron energ y distribution N{E) oc NqE^^ 
with i = 1.5 (Valtaoia et"aLlll992aD. This is shghtlv flatter than 
what i s obtained for' 3 C 273 in iTiirler et af] 12000.) and for 3C 
279 in lLindfors et alj ( l2006h . where values close to s = 2 were 
obtained. In their studies the flare spectra were flattening from 
ttthin = -1-1 to athin = "0.5 as the flare ev olved. This is in 
accord anc e with the original shock model by Marscher & Geail 
(HIH). In lStevens etalTlT994,) the declining slope of -0.28 was 
explained as being due to the radiation being only partially thin 
close to the turnover frequency. This would cause the points to 
deviate from a straight line in the declining part of the spec- 
trum. Also, data at sub-mm frequencies above 230 GHz would 
be required to estimate the accurate value of athin- The rising 
slope of 0.41 is in agr eement with the original shock model 
(Marscher & Gear 1985) and also independently obtained in 
Valtaoia et al. ( 1992a). Therefore we consider our results to be 
in good agreement with the models, considering the scatter in 
our data. 
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6. Conclusions 

In this paper we studied the long-term radio variability of a sam- 
ple of 55 AGNs. We divided the flux density curves into individ- 
ual flares and studied the properties of 159 well-monitored flares 
in at least two frequency bands between 4.8 and 230 GHz. Our 
main results can be summarised as follows: 

1. The flares last on average 2.5 years at 22 and 37 GHz and 
are only slightly longer at lower frequencies. This has im- 
portant consequences when planning multifrequency cam- 
paigns, which usually last only from a couple of days to 
months and therefore are unable to capture a radio flare as 
a whole. 

2. There are no significant diff'erences between the different 
source classes when the durations of the flares are studied al- 
though there is weak evidence for longer flares in the LPQs. 
When Doppler-corrected durations are studied the LPQs still 
have the longest flares but the diff'erence is smaller. Also, 
the BLOs with their shorter flares diff'er from other source 
classes at the lower 4.8-14.5 GHz frequencies. 

3. When studying the duration of the flares against the intrinsic 
Doppler-corrected peak luminosity, we found no correlation 
indicating that the energy release in a flare does not increase 
with duration for these flares. 

4. The flares seem to fo llow the generalised shock model by 
IValtaoia e t al.l (Il992al) quite well even though there is large 
scatter due to incomplete sampling and complicated struc- 
ture of the flares at different frequency bands. 
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Table 1. The sample. For each source the number of flares, the number of datapoints since 1980 in each frequency band, the maximum flux density (S^ax), variability index V and median duration 

of flares at 37 GHz arc listed. Colum 4 indicates if the source is a high-confidence EGRET detection (x) or a possible detection (?). 
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Table 1. Continued. 
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4.0 


o 
O 


1 /I ^ 

14. J 


zz 






IjU Zd\j 
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1 
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~ 


436 


92 


84 
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34 
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X 
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145 


39 
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1 








33 


27 


20 


14 


3.03 
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3.56 
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X 


1 
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909 
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0.62 


2.23 
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X 


6 
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73 
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55 


22 
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' = Values taken at 90 GHz 
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